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Medium-induced broadening and softening of a parton shower 
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The modifications of the angular and transverse momentum distributions of quarks and gluons 
inside a parton shower due to the presence of a medium are studied within an analytical description 
that reduces to the modified leading logarithmic approximation (MLLA) of QCD in the absence of 
medium. 
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The combined measurements of hadrons and photons 
at large transverse momenta performed at the relativis- 
tic heavy-ion collider (RHIC) [H are most naturally in- 
terpreted in terms of the energy degradation of high-pT 
partons in the medium formed in Au-Au or Cu-Cu colli- 
sions prior to their hadronization. A significant amount 
of the fastparton energy is radiated away as soft glu- 
ons 0, [j]: the net effect of the medium is an en- 
hancement of this radiation over the "vacuum" parton 
shower due to the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) parton evolution that leads to jets in 
e+c^ or pp/pp collisions. It is thus expected that the 
characteristics of a "jet" created in an ultrarelativistic 
heavy-ion collision^ differ from those of a jet found in 
collisions of elementary particles. 

Possible medium-induced modifications of jets include 
a broadening of the transverse momentum spectrum 0, S| 
and the distortion of the longitudinal distribution inside 
the jet [§|, increasing jet multiphcities sl. [ol. [lol. 
changes in the jet hadrochemistry [13 L Assessing 
effects is an experimental challenge 1, Kl 
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which would give access to properties of the medium 
probed by the energetic parton. 

The satisfactory description of data on the longitudinal 
and transverse momentum distributions of hadrons in- 
side a "vacuum" jet are successes of perturbative QCD, 
specifically of its modified leading logarithmic approxi- 
mation (MLLA), which takes into account the leading 
double logarithmic terms together with single logarithmic 
contributions , and of next-to-MLLA corrections 

that keep track of energy conservation of parton split- 
tings with increased accuracy [Toj . Therefore, it would 
be desirable to be able to describe in-medium jets with 
a formalism that in the absence of a medium reduces to 
(N)MLLA. In particular, the probabilistic interpretation 
of the latter in terms independent successive branchings 
with angular ordering due to destructive interferences in- 



^ I shall herafter assume that some notion of a jet — whose spec- 
ification might be driven by the attempts at identifying col- 
limated clusters of energy above a fluctuating background in 
data @] — makes sense even in the high-multiplicity environ- 
ment of a nucleus— nucleus collision at high energy. 



side the parton shower makes it highly suited to a Monte- 
Carlo implementation. Since energetic partons created in 
a heavy-ion collision propagate in an expanding medium 
with an anisotropic geometry, a quantitative characteri- 
zation of the resulting partonic shower will require a nu- 
merical approach: the interest of a "medium-modified" 
MLLA is here again obvious. 

While the first numerical implementations of in- 
medium parton showers are appearing 2^ 21 1, it is still 
necessary to have analytical calculations — obviously un- 
der simplified assumptions — of possible properties of the 
modified showers, as references against which Monte- 
Carlo computations can be tested. In this Letter, I 
shall thus present results on the angular and transverse 
momentum distributions of partons inside a medium- 
modified shower. Following the approach proposed in 
Ref. , the influence of the medium is modeled as an en- 
hancement of the singular 1/z parts of the leading-order 
parton splitting functions by a constant factor 1 + /mod 
(this factor is denoted Ns in Refs. l3, 1^): 
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while the other two splitting functions remain unchanged. 
The modeling is admittedly crude, yet has the advan- 
tage of allowing controlled analytical calculations, while 
mimicking the increased radiation of soft gluons of the 
approaches 0, [3, |3| that describe RHIC results on lead- 
ing hadron suppression. Additionally it incorporates the 
conservation of energy at each step of the parton shower 
evolution, which plays possibly an important role in shap- 
ing the global jet characteristics. A more realistic change 
of the splitting functions, which accounts for the coher- 
ent Landau-Pomeranchuk-Migdal radiation in a ther- 
malized weakly coupled deconfined medium can be found 
in Ref. [13] ■ Reference [2^ provides a somewhat equiva- 
lent description in terms of Sudakov form factors. 

Consider a jet with energy E and opening angle Qq. 
For any jet particle with 4-momcntum (/sq, k) one defines 

^ = ln— =ln-, j/ = ln— win — , (3) 

fco X Qa Qo 
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where Qo is the infrared cutoff parameter and k± > Qo 
resp. Q < &o is the transverse momentum resp. angle 
of the radiated parton with respect to the energy flow 
direction ("jet axis"). Finally, let le^ = ln(£'sin9o/(5o), 
corresponding to the maximum transverse momentum of 
a particle in the jet. While approximating sin0 ^ 9 is 
often made for the small opening angles relevant for jets 
in high-energy collisions, in the following I shall retain the 
sine. Additionally, I shall take Qo = Aqcd — 250 MeV, 
which gives a good agreement with data for longitudinal 
hadron spectra, and show results for the values Yoo =5.1 
and Yeo — 6, corresponding to jets with energies typical 
for RHIC (40 GeV) and LHC (100 GeV), respectively. 

The starting point to study the angular distribution in- 
side a jet is the distribution F^^{x, 9, E, 9o) of the energy 
fraction x of hadrons within a subjet with an opening an- 
gle 9 < 9o ■ convolution of the probability 
D4 (u, EQo,uEQ) that the initial parton Aq emit par- 
ton A with the energy fraction u and the virtuality uEQ 
and of the fragmentation function D^Kx/u, uEQ, Qq) for 
the production of hadron h off A with the energy fraction 
x/u and a transverse momentum scale uEQ > Qq: 



Fl{x,e,E,eo) = 

du Dij,u, EOo, uEQ) [^,uEQ, Qo) ■ (4) 
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The number of hadrons h inside the opening angle 9 is 
given by the integral of this distribution over x [2J] : 

Nl{Q,E,Q^)^ 
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du uD^^j^u, EQo, uEQ) N^{uEQ, Qo), (5) 



where N^{uEQ, Qo) denotes the multiplicity of hadrons 
inside the "subjet" initiated by parton A with the vir- 
tuality uEQ. The physical picture of Eqs. (|3][5]) is that, 
as a consequence of the angular ordering of the splitting 
processes, the hadrons that end up in the subjet origi- 
nate from quarks and gluons emitted in the shower when 
its typical fcj_ scale was uEQ, after an evolution from the 
initial scale down to that intermediate scale. 

Since the purpose in this Letter is to establish the mod- 
ifications of jets induced by a medium compared to their 
"vacuum" MLLA properties, for the hadronic fragmen- 
tation function relevant for vacuum jets I shall use 
the corresponding MLLA (partonic) function D^^^^, in 
the limit Qo Aqcd- The latter is in the case of a 
gluon jet related to the MLLA "limiting spectrum" (l7j 
through L>i™(£,y) = xDf^^^{£,y). Note that exact en- 
ergy conservation is spoiled in the approximations lead- 
ing to D^^^^ and thus Z)'""; yet the latter still represents 
a good baseline at sufhciently small x. Consequently, the 
hadron multiplicity inside a gluon jet is given by 
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FIG. 1: Angular dependence of multiplicities inside gluon jets 
with E = 100 GeV (Ygo = 6) in MLLA and for f^cd = 1, and 
with £ ~ 40 GeV = 5.1) in MLLA and for /^^d = 0.8. 

and that inside a quark jet by (25j 




Ki^e), (7) 



where Yq = ln(£;sin9/Qo), a = {IIN^ + 2Nf/N^)/3, 
(3o = [IIN^ - 2iV/)/3 and B = a//3o, while /C'^ is, follow- 
ing the local parton-hadron duality (LPHD) hypothesis, 
a proportionality constant relating partonic and hadronic 
properties. The actual value of K/'' will be unimportant 
in the following, since I shall assume that the same value 
should hold for vacuum and in-medium jets. 

The longitudinal ln(l/a;) distribution of partons inside 
a (gluon) jet within the adopted modeling of medium 
effects was reported in Ref. j9| , and amounts to a strong 
distortion of the MLLA hump-backed plateau. This leads 
to increased (sub)jet multiplicities, /So in Eqs. being 
replaced by /3o/(l + /mod) and a by a H- 4iVc/mod- 

Regarding , for vacuum jets it is given by the 
usual expression following from the Altarelli-Parisi par- 
ton splitting functions. Replacing the latter with the 
modified functions (IT][2]) and performing a Mellin trans- 
formation as well as an expansion of around u = 1 in 
Eq. ([5]), one obtains readily the multiplicity A^^,^ which 
is plotted vs. the opening angle 9 in Fig. [1] for both 
MLLA and medium-distorted jets. The value /mcd = 0.8, 
which allows to reproduce a nuclear modification factor 
Raa ~ 0.2 [§], was adopted for the RHIC energy, while 
/med = 1 was choscu for the higher (LHC) energy. As in 
Ref. 'iil, the LPHD constant was taken to be J& = 0.2. 
Figure [T] shows that the gradient diV^ /d9 is larger over 
a broader angular range for the medium-modified than 



^ The detailed calculations will be reported elsewhere. 
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for the MLLA gluon jet. This behavior, which is also ob- 
served on quarks jets (not shown here), is the in-medium 
angular broadening of the parton shower which was ob- 
served in Q-Pythia [uj. Note that in Fig. [T] the low-9 
region has been suppressed, since even the MLLA results 
are not reliable for small angles (see the remark on low 
k±_ values below). 

The hadron multiplicities represented in Fig. [1] in- 
clude all hadrons down to the infrared cutoff Qq. Given 
the high-multiplicity environment of heavy-ion collisions, 
experimental studies will most probably rather investi- 
gate the modifications of jet properties above some pr 
(with respect to the beam axis) cutoff which might be 
quite large. In that case, the angular broadening is sig- 
nificantly less marked.^ 

Let me turn to the distribution of transverse momenta 
k± in a jet. Deriving Eq. ^ with respect to ln0 — or 
equivalently, at fixed x, In k± — yields the double differen- 
tial single-particle inclusive distribution of hadrons inside 
a jet with opening angle 80 and energy E [l7l |: 



da; d In fc I 



dxdlne dine 



Fl{x,e,E,eo). (8) 



Within a small-x approximation, one can expand 
around u = 1 in Eq. This allows one, after some 

algebra [2^, to rewrite in terms of the average color 
current {C)ao inside a jet initiated by parton Aq and of 
the distribution in ^ = ln(l/a;) of hadrons inside a gluon 
jet, 'Dg{£,y) = xDg{x, EQ,Qo). Inserting the result in 
Eq. dH) yields 





d 




dy 



(9) 



This expression was computed for vacuum jets within 
MLLA [2^, using the limiting spectrum for Vg, and 
later including some next-to-MLLA corrections, which 
yield improved agreement with measurements by the 
CDF Collaboration for k± > 1 GeV/c 0. Given the 
exploratory approach pursued in this Letter, I shall only 
investigate the influence of the medium modifications in 
Eqs. ((TJH]) on the MLLA fcj^-distributions. The average 
color currents {C)g^q following from the modified splitting 
fimctions can be computed^ and inserted in Eq. (l9|) to- 
gether with the distorted longitudinal spectrum [§]. For 
values of £ where the assumptions behind the MLLA com- 
putation are fulfilled (see discussion in Ref. [1^), the 
resulting double differential single-particle distribution 
is, with respect to the vacuum one, depleted at low fc^ 
{y ^ 1.5) and significantly larger at higher k± 27[. The 
growth at high k±_ obviously reflects the angular broaden- 
ing already observed in Fig.[Tl while the low-fc_L behavior 
should probably be considered with care, since it corre- 
sponds to a region where MLLA is not to be trusted. 
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FIG. 2: In fcx-distribution inside gluon jets for the same cases 
as in Fig. [T] 

Integrating Eq. ([9]) over i yields the transverse momen- 
tum distribution inside a jet: 



dA^ 
d In /s I 



d^ d In A: I 



(10) 



The latter is displayed in Fig. [2] for vacuum and medium- 
distorted gluon jets at RHIC and LHC energies. While 
the values below y = ln(A:^/(5o) = 1-5, should not be 
taken too literally — the divergence at y ^ 0, which re- 
flects that of the running QCD coupling constant as {k± ) 
when k± — *■ Aqcd, hints at the breakdown of the per- 
turbative regime — , the overall trend of medium effects 
seems to be a push of the distribution towards larger 
transverse momenta: this is a transverse momentum 
broadening of the jets. Note that the "usual" source of 
transverse momentum broadening (which was however 
not seen in the first Jewel simulations |2j]) is the trans- 
fer of momentum to the jet in the interactions between 
the latter and the medium 0: this is not taken into 
account in the model. The phenomenon reported here 
is different, since no momentum is injected into the jet 
from the outside. The eventual broadening of the jet will 
result from both effects, which do not necessarily add 
up constructively. Investigating that issue is beyond the 
scope of the present study 

Since the effect of the medium is to redistribute par- 
tons from high towards low x- values 0], it was not to- 
tally obvious that this would at the same time involve a 
wider k± distribution as seen in Fig. [21 Both modifica- 
tions can simultaneously happen only if the presence of 
the medium makes the transverse momentum spectrum 
steeper. This is illustrated in Fig. [31 which indeed shows a 
medium-induced depletion of the yield at high fcj^ , which 
was barely visible on Fig. [21 together with the enhance- 
ment at low k±. All in all, the k±_ spectrum is softened, 
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well as a softening of the k± spectrum. The multiplicity 
over a lower k± cutoff, which is more easily accessible ex- 
perimentally, shows an enhancement in medium-modified 
jets up to large cutoff values. 

Strictly speaking, the results obtained here hold at 
the parton level. Using the LPHD hypothesis, which for 
(N)MLLA jets leads to a good agreement with measure- 
ments of ln(l/x) and k± spectra, they would translate 
into characteristics at the hadronic level. The hadroniza- 
tion process might however possibly totally reshuffle the 
distributions inside the jet [_2C/|, in which case the predic- 
tions could only be checked in Monte-Carlo simulations. 



I wish to thank Redamy Perez-Ramos for enlightening 
discussions. 



FIG. 3: Transverse momentum distribution inside gluon jets 
for the same cases as in Fig. [1] 




FIG. 4: Ratio of the multiplicities above a given k± cutoff 
inside medium-distorted and vacuum jets, for gluon and quark 
jets at typical LHC (Veo = 6) and RHIC CVeo — 5.1) energies. 

which was also seen in Q-Pythia simulations 

Integrating the k± spectrum above a given A;™* yields 
the multiplicity above that cutoff, which might be eas- 
ier to assess experimentally. One can then form some 
"medium modification factor" by dividing the multiplic- 
ity in a medium-distorted jet by that inside a vacuum 
jet. This ratio is shown for both gluon and quark jets in 
Fig. H it is larger than 1 up to fc™* ~ 10 - 15 GeV/c, 
and is systematically larger for quark than for gluon jets. 

Using an analytical formalism, I have computed the an- 
gular and transverse momentum distributions of partons 
inside a medium-distorted parton shower. The effect of 
the medium is, as anticipated, a broadening of the distri- 
butions with respect to those of "vacuum" MLLA jets, as 
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